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ABSTRACT 



Context. Magnetohydrodynamic thermal modes may play an important role in the formation, plasma condensation, and evolution of 
solar prominences. Unstable thermal modes due to unbalance between radiative losses and heating can lead to rapid plasma cooling 
and condensation. An accurate description of the radiative loss function is therefore crucial for this process. 

Aims. We study the stability of thermal modes in unbounded and uniform plasmas with properties akin to those in solar prominences. 
Effects due to partial ionization are taken into account. Three different parametrizations of the radiative loss function are used. 
Methods. By means of a normal mode analysis, we investigate linear nonadiabatic perturbations superimposed on the equilibrium 
state. We find an approximate instability criterion for thermal modes, while the exact linear growth rate is obtained by numerically 
solving the general dispersion relation. The stability of thermal disturbances is compared for the three different loss functions consid- 
ered. 

Results. Using up-to-date computations of radiative losses derived from the CHIANTI atomic database, we find that thermal modes 
may be unstable in prominences for lower temperatures than those predicted with previously existing loss functions. Thermal insta- 
bility can take place for temperatures as low as 15,000 K, approximately. The obtained linear growth rates indicate that this instability 
might have an important impact on the dynamics and evolution of cool prominence condensations. 

Key words. Instabilities — Sun: filaments, prominences — Sun: corona — Sun: atmosphere — Magnetohydrodynamic s (MHD) 



1. Introduction 

Thermal or condensational modes have been extensiv e ly in- 
vestigated in magnetiz ed plasmas (e.g., Parker] 1 1953t iFieldl 
1965; lHevvaertsl[l974l) . As explained by Parkel dl953h . ther- 
mal instability can happen in a diffuse medium due to un- 
balance between temperature-independent energy gains , i.e., 
heatin g, and temperature-dependent radiative losses. iParketl 
d!953l) arrived at the qualitative criterion that instability can 
be present when radiative losses decrease as the tempera- 
ture increases. This mechanism is important in the context 
of prominences since it allows the formation of cool plasma 
condensations in a medium of high temperature. Therefore, 
unstable thermal modes may play an important role in the 
formation of solar promi n ences and in the evolution of the 
prominence plasma. IFieldl d 19651) and iHevvaerts ( 19741) in ves- 
tigated the phenomenon in more detail than Parker (119531) and 
derived more accurate instability criteria. Subsequent papers 
that investigate both linear and nonlinear ther mal instabilities , 
mainly in the field of prominences, are, e.g., Hildnet (Il974l): 
lOranetalJ dl982l): iDahlburg & Mariskal d!988l): iKarpen et alJ 
1989 1) : ICargill & Hoodl dl989h: IVan der Linden & Goossensl 



1991 



): ICarbonell et all d2004l) : 



Soler et al.1 (1201 lh among oth- 



ers. 



As shown in the instability criterion derived bv lFieldl d 19651) . 
an accurate description of the radiative loss function is crucial 
to ascertain the stability of thermal modes. However, the deter- 
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mination of the radiative loss function in prominence plasmas 
depending on the values of temperature and density is a difficult 
work that requires complicated numerical solutions of the radia- 
tive transfer equations for nonlocal thermodynamic equilibrium. 
Alternatively, several semi-empirical parametrizations of the ra- 
diative loss function for prominence and coronal conditions are 
available in the literature. These parametrizations enable us to 
incorporate radiative losses in a consistent way in the theoretical 
models of prominence plasmas without the need of solving the 
radiative transfer problem. 

One of t he firs t parametrizations of this kind can be found 
in iHildneti d 19741) . who performed a piecewise fit as func- 
tion of the temperature of the computations of radiative losses 
available by then. Subsequent authors have pro posed differ- 
ent parametrizations that update Hildner's fit (e.g. jRosner et al.l 
11978b iMilne et all Il979h . although Hildner's function is still 
used in some works nowadays. Another function frequently 
used in the lite rature is the so-called Klim chuk-Raymond func- 
tion (see, e.g.. iKlimchuk & CargilJ l200lh that is a better rep- 
resentation of the radiative losses for prominence-corona tran- 
sition region and coronal temperatures. The shape of the loss 
function depends on the completeness of the atomic model 
used for the calculation, on the atomic processes included, on 
the ionization equilibrium, and element abundance assumed. 
More recent loss functions that incorporate accurat e atomic data 
inform atio n are, e.g., the loss f unction used by iParenti et al.l 
d2006l) an dlParenti & Viall d2007l) computed from the CHIANTI 
database (Dere et alJll997tlLandi et alJl2012l) . and the loss func- 
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lion computed by Schure et all (1201 Ol) using the SPEX package 
(IKaanstra et al.lll996l) . 

The purpose of this paper is to investigate the stability prop- 
erties of thermal modes in cool prominence plasmas. We com- 
pare the results using up-to-date computations of radiative losses 
derived from the CHIANTI v7 atomic database with those ob- 
tained assuming two of the most used loss functions existing in 
the literature. These radiative losses are obtained assuming an 
optically thin plasma, while the core of prominence cannot com- 
pletely satisfy this condition. We will discuss our assumption 
later on. This paper is organized as follows. Section [2] contains 
a description of the equilibrium and the basic equations. The 
instability criterion for thermal modes is derived in Section [3] 
Parametric studies of the linear growth rate are done in Section[4] 
Finally, Section [3] contains the summary and discussion of the 
results. 

2. Basic equations 

2.1. Equilibrium 

Our equilibrium configuration is a uniform plasma of infinite ex- 
tend. We assume a partially ionized hydrogen plasma composed 
of ions, electrons, and neutrals. We use Cartesian coordinates 
and all quantities are expressed in MKS units thorough this pa- 
per. The magnetic field, B, is uniform and orientated along the 
z-direction, namely B = Be z , with B constant. We denote by p, T, 
and p the equilibrium mass density, temperature, and gas pres- 
sure, respectively. The set of basic nonadiabatic MHD equations 
governing the plas ma dynamics in the single-fluid approxima- 
tion are (see, e.g.,[Braginskii 1 965I) 

Dp 

Df 



-pV • v, 



(1) 



Dv 1 
p — = -Vp + -(VxB)xB, 
Df p 

3B 

— = Vx(vxB)-Vx(^VxB) 
dt 



Vx 



I B 2 



[(V x B) x B] x B 



^ = IP^ + iy - m . iKm - P L(p,T)]. 
Df p Df 



P = 



pRT 
P 



(2) 



(3) 



(4) 



(5) 



where j=| = Jp + v • V is the material derivative for time vari- 
ations following the plasma motion, v is the plasma velocity, p 
is the magnetic permittivity, y is the adiabatic index, k is the 
thermal conductivity tensor, L (p, T) is the heat-loss function, rj 
and rjc are the coefficients of Ohm's and Cowling's diffusion, 
respectively, R is the ideal gas constant, and p = is the 
mean atomic weight, with £ the relative fraction of ions. This 
parameter ranges from £ = 1 in a fully ionized plasma and 
£ = in a neutral gas. In Equations (U)-©, the effects of grav- 
ity and viscosity have been omitted. Equation (O is the induc- 
tion equation. In Equation (f3]) we have neglected some minor 
terms which are several orders of magnitude smaller than Ohm's 
and Cowling's terms in partially ionized prominence plasmas 
(see an expr ession of the complete single-fluid induction equa- 
tion in, e.g.. iForteza et al.l r2007t IZaqarash vili et al.ll201 lb . The 
neglected terms are Hall's ter m, the diamagnetic current term, 
and Biermann's battery term. ISoler et all (l2009bh showed that 
the effect of Hall's term on the waves is negligible in promi- 
nence conditions. The diamagnetic current term and Biermann's 



battery term are only relevant when large pressure gradients are 
present, a situation more representative of stellar interiors. The 
omission of these terms in the present work is therefore justified. 
In a partially ionized plasma, Cowling's diffusion represents an 
enhanced magnetic diffusion due to ion-neutral collisions. The 
expressions for rj and rjc are 



rj = 3.7 x 10 



_ 6 m e In A c „_ 



3/2 



pe^- 

£ 2 (l-£) 2 

pa n 



(6) 



(7) 



where m e is the electron mass, e i s the electro n charge, In Ac 
is Coulomb's logarithm (see, e.g.. lPries3 [T982). and a n is the 
neutral friction coefficient given by 



a n = -6(1 -£) — • 



(8) 



with m p the proton mass and kg Boltzmann's constant. In the 
fully ionized case, rjc — r\ and the third term on the right-hand 
side of Equation (0 is absent. 

We denote by K\\ and k± the parallel and perpendicular scalar 
components of the thermal conductivity tensor with respect to 
the magnetic field direction, which can be expressed as 



*|| = K e + K B , 



K ± =Ki + K n 



(9) 



with K e , k\, and a:,, the conductivities by electrons, ions, and neu- 
trals, respectively. In a fully ionized medium, K\\ is governed by 
electrons, whereas k ± is caused by ions. In a partially ionized 
plasma, the contribution of neutrals, A: n , has to be added to both 
scalar conductivities because thermal conduction by neutrals is 
isotropic. Since k\ <k K n , the conductivity by io ns can be ne- 
glected in k ± . We use (see, e.g JSoler et al.| [2010) 



K e = Kr n 6r 5/2 , 



* n =2.24x 1(T 2 (1 



.1/2 



(10) 



Regarding the heat-loss function, L (p, T), we consider the 
following expression, 



L(p,T)= PX *T a -h, 



(ID 



where x* ar, d a are functions of the temperature, and h is an 
arbitrary heating function. In the equilibrium we assume that 
radiative losses balance heating, so that h is defined to satisfy 
L (p, T) = in the equilibrium state. Several parametrizations of 
X* and a for prominence conditions are available in the litera- 
ture. In this work, we use two of the most used loss f unctions, 
name ly the well-known paramet rizations by Hildner (IHi ldner 
11974 and Klimchuk-Raymond dKlimchuk & Cargillll200ll) . In 
addition, we use a parametrization of the radiative l oss function 
comp uted from the CHIANTI v7 ato mic database dLandi et aT] 
120121) assuming coronal abundances (Mazzotta et al. 119981) . ion- 
ization equilibrium, and a constant pressure of 6.64 mPa. The pa- 
rameters x* an d <* corresponding to these three different fits are 
given in Table [1] Figure Q]displays the loss rate per unit volume 
as function of the temperature for the three fits. The main differ- 
ences in the three curves are found at low (T < 3 x 10 4 K) and 
high (T > 10 6 K) temperatures. For the region of our interest, 
that is the low temperature range, the peak at about 1.5 x 10 4 K 
in the CHIANTI calculation is due mainly to the H and He 
emissions. Besides this difference, it should be noted that the 
CHIANTI database is still incomplete at these temperatures, so 
that further increase of the loss function is expected for more 
complete calculations. 
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Table 1. Parameters of the loss function for the considered fits. 
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Fig. 1. Fit of the energy loss rate per unit volume as func- 
tion of the temperature computed from the CHIANTI atomic 
database assuming coronal abundances, ionization equilibrium, 
and a constant pressure of 6.64 mPa (solid line) compared to 
Hildner's (dotted line) and Klimchuk-Raymond's fits (dashed 
line). The parameters x* an d a of the three fits are given in 
Tabled 

2.2. Dispersion relation for linear perturbations 

We take the plasma initially at rest and superimpose linear per- 
turbations on the equilibrium state. Equations (H)-© are lin- 
earized. We write perturbations proportional to exp (/k • r - iojf), 
where r = (x,y,z) is the position vector, k = (k x ,ky,k^j is 
the wavenumber vector, and oj is the frequency. For simplic- 
ity, we choose the reference frame so that we can set k y = 
and consider wave propagation in the xz-plane only. We focus 
our study on nonadiabatic magnetoacoustic and thermal modes. 
Alfven modes are not discussed in the present investigation. We 
combine the linearized Equations (Q~|)-(0 and, after some al- 
gebraic manipulations, we arrive at the dispersion relation for 
nonadiaba tic magnetoacoustic and thermal mod e s (see details 
in, e.g., ICarbonell et ail 12004 ISoler et alj I20T0I: IBarcelo et all 
which can be written in a compact form as follows, 

w 4 _ ( r 2 + A 2) ^2^2 + k 4 A 2 ^2 _ y 2 ) + y 2 cofj 2 Q j = ^ (u) 

where k 2 — k 2 + k 2 is the square of the wavevector modu- 
lus, cos# = kjk is the cosine of the angle between k and B, 
A 2 = yp/p is th e square of the nonadiabatic sound speed (see 
ISoleretalJl2008h . with y the effective nonadiabatic index de- 
fined as 



r = 



with 



(r -!)[(■ 



k ± sin 2 6 + K\\ cos 2 6*) k 2 + o>j - o» p J - 



iyo> 



(y - 1) [(k ± sin 2 9 + k\\ cos 2 k 2 + w T J - ioj 



T T 

K\\ = —K\\, Kj_ = —K±, 



to = P -(— 

P P \dp) T 



pL{— 

p\dT 



(13) 



(14) 



(15) 



In addition, Y 2 is the square of the modified Alfven speed (see 
ISoler et al. 2009a) defined as 
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(16) 



Notes. The fit derived from computations based on th e CHIANTI 
atomic database assuming coronal abundances was used in Paren tiet al.l 
i2006h. Hildn er's and Klimchuk-Raymond's fits are adapted from 
Hildner ( 1971) and lKlimchuk & CargililfcOOlh . respectively. Quantities 
are expressed in MKS units. 



with v 2 = B 2 lnp the square of the ideal Alfven speed. Note 
that both A 2 and T 2 are functions of to and k. In the absence of 
nonadiabatic effects, x|| = ic ± = coj = io p — and y — y, so 
that A 2 becomes the square of the adiabatic sound speed, v 2 = 
yplp. In the absence of magnetic diffusion, rjc — r\ — and T 2 
becomes the square of the ideal Alfven speed, v 2 K . In such a case, 
Equation Sl2\ reverts the the well-known dispersion relation for 
ideal, adiabatic magnetoacousti c waves in a plasma of infinite 
extend (see, e.g., Lighthill 196XJ)- 

For positive and real k and 8, Equation (fT2l i has five solutions 
of oj. In general, u> is complex, namely co = oya+ioji, with <wr and 
oj\ the real and imaginary parts of co, respectivel y. The solutions 
of Equ ation (TTZll were di s cussed in detail by Carbonell et al.l 
(120041) and IBarcelo et ail (1201 lb . Of the five solutions of oj, 
two complex conjugate solutions correspond to damped slow 
modes and other two complex conjugate solutions correspond 
to damped fast modes. The remaining solution is purely imag- 
inary, i.e., wr = 0, and corresponds to the thermal mode. The 
thermal mode is the subject of our investigation. 



3. Approximate instability criterion 

We perform a first-order expansion of Equation ( TTZl i for a low-/? 
plasma, where (3 is the ratio of the gas pressure to the magnetic 
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Fig. 2. Approximate thermal mode growth rate vs. tempera- 
ture in the absence of thermal conduction and for Hildner's 
(dotted), Klimchuk-Raymond's (dashed), and CHIANTI-based 
(solid) loss functions. The shaded area denotes the region of in- 
stability at low temperatures obtained with the CHIANTI-based 
loss function that is not present for the other parametrizations. 



pressure. We obtain two different approximate dispersion rela- 
tions, namely 



aV 

oj 2 - k 2 T 2 - — — ^-k 2 sin 2 0*0, 
F 2 - A 2 

for fast modes, and 



■k 2 A 2 - 



A 



2 v 2 



t 2 -a 2 



k 2 sin 2 « 0, 



(17) 



(18) 



for slow and thermal modes. Fast modes are weakly affected by 
nonadiabatic mechanisms and are not investigated further in the 
present work. We focus on Equation ( TT8l and explore in more 
detail the approximation for thermal modes. Note that in the low- 
f3 regime the propagation of slow and thermal modes is almost 
parallel to the magnetic field. By assuming F 2 » A 2 , i.e, the 
Alfven speed is much larger than the sound speed as typical for 
prominence conditions, Equation (fT8l simplifies to 



■k 2 A 2 



1 



v 2 



0. 



(19) 



Now we use the definitions of A 2 and T 2 and expand 
Equation (fT9l as a third order polynomial in u>. Since the thermal 
mode is a purely imaginary solution we write a> = wr+is, where 
s is the thermal mode growth rate, and set wr = 0. Because 
of the temporal dependence exp (—icoi), thermal mode perturba- 
tions are proportional to exp(sf), meaning that for s > pertur- 
bations grow in time. This behavior corresponds to instability. 
On the contrary, for s < thermal disturbances are damped. To 
obtain an approximation for s, we neglect terms with <3(s 2 ) in 
the polynomial expansion of Equation ( fT9] l. The neglected terms 
are related to the slow mode and produce minor corrections to 
the thermal mode growth rate. After some algebraic manipula- 
tions, we obtain the approximate growth rate, namely 



Uk± sin 2 8 + k\\ cos 2 6j k 2 + lo-y - w p J v 2 cos 2 8 
[(k± sin 2 8 + K\\ cos 2 dj k 2 + cjj - w p J j] C + y ' A ^ 



(20) 



In the absence of magnetic diffusion, rjc 
becomes 



and Equation ( 1201 



r ■ 



k ± sin 8 + ki I cos 



8)k 2 



(21) 



Equation (f2Tb is similar to th e e xpressions found by 
IVan der Linden & Goossensl (1 19911) and ISoler etalj d2011l) for 
thermal con tinuum modes, and by ICarbonell et al. (2009J) and 
ISolerl J20T0I) for the imaginary part of the frequency of propa- 
gating thermal waves in a flowing medium. In the case without 
thermal conduction, i.e., K\\ — k± = 0, the approximate growth 
rate is 



y-i 

y 



(22) 



Equation (l22l is independent of k, 8, and the magnetic field 
strength and orientation. 

Equation (120b provides us with the instability criterion. 
By taking into account the definitions of % a>j and <d p 
(Equations ( fl4b and (fT3l)). we find that the combination of pa- 
rameters to have s > in Equation ( f20b must satisfy the condi- 
tion 



y - 1 v 2 A cos 2 ( 



— < (k ± sin 2 8 + k\\ cos 2 6*) k 2 



dL 
BT 



P 



T \ dp) T 



(23) 



To the best of our knowledge, Equation d23l is the first in- 
stability criterion for thermal modes that accounts for the ef- 
fect of Cowling's diffusion. We deduce from Equation ( 1231 that 
Cowling's diffusion has a stabilizing role since Cowling's dif- 
fusion incorporates a lower bound in the instability criterion. In 
order to compare with previous instability criteria derived in the 
literature, we set 77c = and Equation (l23l simplifies to 



(k x sin 2 8 + K\\ cos 2 8^k 2 +p y^j 



81 
T 



dp) T < °' 



(24) 



Equat i on (|24l) agrees with the instability criterion provided by 
iFieldl (1 19651) in his Equation (25a) if the term accounting for 
perpendicular thermal conduction in added to Field's expression. 
For simplicity, we take the case r]c — and use the expression of 
the heat-loss function L (Equation (fTTT i) to rewrite Equation (124-b 
in terms of parameters x" an d a - Then, instability is present for 
values of a satisfying 



a < 



\K X sin 2 8 + K\\ cos 2 6>) k 2 
1 p 2 X *T a - x 



(25) 



In the absence of thermal conduction (k\\ — k ± - 0) or for long 
wavelengths (k — > 0), the second term on the right-hand side 
of Equation d25l l vanishes and the instability criterion reduces to 
a < 1 . Taking into account the values of the parameter a given in 
TableQ] the lowest thermally unstable temperatures according to 
this criterion are T ~ 8xl0 4 Kin Hildner's fit, T « 9.33 x 10 4 K 
in Klimchuk-Raymond's fit, and T » 1.58 x 10 4 K in the 
CHIANTI-based fit. Importantly, we find that the threshold tem- 
perature for the thermal instability is substantially reduced using 
the CHIANTI-based radiative losses in comparison to Hildner's 
and Klimchuk-Raymond's functions. 

Let us do a simple computation in the case without ther- 
mal conduction (k\\ — k± = 0) and without Cowling's diffu- 
sion (i]c = 0). In this case the approximate growth rate is given 
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by Equation (l22l >. Figure [2] displays the growth rate computed 
from Equation (1221 as a function of the temperature. A constant 
gas pressure of 6.64 mPa is assumed and the density is com- 
puted accordingly. As predicted by the instability criterion, there 
is a region of instability at low temperatures obtained with the 
CHIANTI-based loss function that is not present for the other 
parametrizations (see the shaded area in Figure 0. This region 
of instability is present for temperatures between 1.58 x 10 4 K 
and 3.16xl0 4 K. At these temperat ures the prominence plasma is 
only partially ionized (see, e.g.. lGouttebroze &~L abrosse 2009) 
and the roles of thermal conduction by neutrals and Cowling's 
diffusion may be relevant. These effects are investigated numer- 
ically in the Section|4] 

Additionally, in Figure [2] we notice the abrupt jumps of the 
growth rate at the boundaries of the regions where different val- 
ues of x* an d a are used. The reason of these jumps is that, al- 
though the cooling function L (p, T) is continuous, its the deriva- 
tives with respect to density and temperature are discontinuous 
where the values of parameters y * and a change. Additiona l 
co mments on this issu e are given in [Van der Linden et al. I (fl99l 
and lSoler etal] d2011l) . 

4. Numerical results 

Here we compute the thermal mode growth rate by solving 
the general dispersion relation (Equation ( fT2b ) by standard nu- 
merical methods. We focus on the region of instability at low 
temperatures obtained with the CHIANTI-based loss function 
(shaded area in Figure|2]). In the following computations we use 
B = 10 G, p = 6.64 mPa, and 9 = n/4. These values of the equi- 
librium magnetic field stren gth and pressure are typical promi- 
nence parameters (see, e.g. jLabrosse et ai] i2010). We compute 
the growth rate, s, as a function of the perturbation wavelength, 
A = 2n/k. Wavelengths typically observed i n prominences are 
rough ly between 10 3 km and 10 5 km (see lOliver & Ballesterl 
2002). This range of wavelengths correspond to disturbances 
usually detected from Doppler time series. We do not claim that 
all the observed disturbances are thermal modes. It is very diffi- 
cult to determine the nature of the waves in the absence of addi- 
tional information as, e.g., the velocity polarization or the mag- 
netic and density perturbations. In particular, a clear distinction 
between slow and thermal modes may be very c omplicate (see 
a discussion on this issue in Carb onell et alJl2009l). In this paper 
we us e the observed wavelengths given by lOliver & Ballesterl 
(2002) as the most probable range of wavelengths for thermal 
modes. 

4. 1 . Fully ionized plasma 

First, we investigate the fully ionized case and set £ = 1. 
Therefore, the effect of thermal conduction by neutrals is ab- 
sent and Cowling's diffusion becomes classical Ohm's diffusion, 
i.e., Tjc - rj. FigureOa) displays s versus A. As we focus on the 
unstable behavior of the thermal mode, only positive values of s 
are displayed in Figure|3la). We consider three different temper- 
atures within the region of instability denoted by a shaded area 
in Figure [2] We compare the numerically computed growth rate 
(solid lines) with the approximation of Equation (f2TJ in the ab- 
sence of diffusion, i.e., 77 = (dashed lines). The approximate 
result is in good agreement with the actual growth rate. Different 
test computations with and without Ohm's magnetic diffusion 
(not displayed here for simplicity) indicate that magnetic diffu- 
sion has almost no impact on the value of the growth rate. This 
means that we can safely neglect the effect of diffusion and use 



the approximate growth rate given in Equation (f2TT > for the case 
77 = 0. 

Figure[2a) also shows that the thermal mode is stabilized for 
short wavelengths. The stabilization is due to thermal conduc- 
tion by electrons. Thermal conduction becomes important as the 
wavelength decreases, so that the growth rate is reduced and the 
instability is suppressed for short enough wavelengths. For long 
wavelengths the growth rate saturates and becomes independent 
of A. 

Regarding the temperature, we find that the growth rate de- 
creases as the temperature increases within the region of instabil- 
ity. For T = 16, 000 K the maximum growth rate is s « 0.6 s , 
while for T = 30, 000 K the maximum growth rate decreases 
to s « 0.05 s _1 . These values of the linear growth rate indi- 
cate that the thermal instability operates in short timescales and 
suggest that the effect of the instability in prominences may be 
observable on the dynamics and evolution of cool plasma con- 
densations. Nonlinear studies beyond the present linear analysis 
are needed in order to asses the actual impact of the instability 
on the condensation dynamics. 

4.2. Partially ionized plasma 

Here we incorporate the effects of partial ionization, namely 
thermal conduction by neutrals and Cowling's diffusion. In 
these computations we fix the temperature to T — 16,000 K. 
Figure |3lb) displays s versus A for different values of First of 
all, we obtain that when £ decreases, the critical wavelength for 
stabilization increases due to thermal conduction by neutrals. At 
cool prominence temperatures, thermal conduction by neutrals 
is more efficient than conduction by electrons. Thus, the critical 
wavelength for stabilization is now determined by the conduc- 
tivity of neutrals. We also see in Figurejjb) that the growth rate 
increases when £ is reduced. The growth rate for £ = 0.1 is 
about four times larger, approximately, than that for £ = 1 . We 
can explain these result by taking into account that a constant gas 
pressure is assumed in the computations. So, when £ decreases, 
the effective plasma density grows because of the increase of the 
amount of neutrals. The increase of the growth rate is a conse- 
quence of the increase of the effective density. 

Also, we compare the numerical results (solid lines in 
Figure [Hb)) with the approximate growth rates (dashed lines) 
given by Equation ( f2TT >. Note that Equation (l2TT i misses the effect 
of Cowling's diffusion. Nevertheless a reasonably good agree- 
ment between both results is obtained. The differences get larger 
when £ becomes small. As happens for Ohm's diffusion, we 
find that Cowling's diffusion have little influence on the ther- 
mal mode growth rate. Although Cowling's diffusion imposes a 
lower bound for the instability criterion (see Equation (l23l). its 
influence on the growth rate is of almost no relevance for realis- 
tic values of rjc- 

5. Discussion 

In this paper we have investigated the stability of thermal 
modes in partially ionized prominence plasmas in the single- 
fluid approximation. We have restricted ourselves to the lin- 
ear phase and have derived an instability criterion that takes 
into account the effects of thermal conduction by electrons and 
neutrals, and Cowling's diffusio n. We have a pplied the insta- 
bility criterion using Hildner's (|Hildnerl Q974I) and Klimchuk- 
Raymond's dKlimchuk & Car gill 2001) loss functions, which 
are frequently used in the literature, along with a n ew loss func- 
tion derived from the CHIANTI atomic database dParenti et al.l 
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Fig. 3. (a) Numerically computed growth rate (solid) and approximate value obtained from Equation (fJTJ (dashed) vs. wavelength 
for the CHIANTI-based loss function and different values of the temperature in the case of a fully ionized plasma, i.e., £ = 1. (b) 
Same as panel (a) but for T = 16, 000 K and different values of 



2006] lParenti& Vialll2007h . Results using both Hildner's and 
Klimchuk-Raymond's loss functions predict the presence of 
thermal instability for temperatures higher than 10 5 K, approx- 
imately. However, the threshold temperature is significantly re- 
duced for about an order of magnitude when the more up-to-date 
function based on the CHIANTI database is used. In particular, 
instability can happen at temperatures as low as 15,000 K, ap- 
proximately. Effects due to partial ionization, specially thermal 
conduction by neutrals, are relevant at these low temperatures. 

Focusing on the region of instability at low temperatures 
obtained with the new CHIANTI-based loss function, we have 
performed a parametric study of the linear growth rate by nu- 
merically solving the general dispersion relation. For constant 
gas pressure, we find that the growth rate decreases as the tem- 
perature increases. In addition, the growth rate increases as the 
amount of neutrals gets larger. We also find that thermal conduc- 
tion reduces the growth rate for short wavelengths, with conduc- 
tion by neutrals being more efficient than conduction by elec- 
trons in partially ionized plasmas. However, neither Ohm's nor 
Cowling's diffusion have an important influence on the growth 
rate. 

A few remarks should be made about the applicability of our 
results. We must note that assumption of optically thin plasma 
might no be valid at low prominence temperatures. Due to fi- 
nite optical thickness, the actual radiative losses of the plasma 
would be reduced in comparison to the optically thin case, so 
that t he thermal mode gr owth rate would decrease consistently 
(see lCarbonell et aI1l2006l) . Some attempts to incorporate the ef- 
fect of finite optica l thickness in Hildner's parametr i zatio n can 
be found in, e.g., iRosner et al.1 d!978l) : iMilne et alj d 19791) . To 
our knowledge, the effect of finite optical thickness has not been 
incorporated in more up-to-date loss functions. We recall that 
the CHIANTI database is still incomplete at temperatures of the 
order ofT~ 10 4 K. Further increase of the loss function at low 
temperatures is expected for more complete calculations. Thus, 
these two effects, namely the decrease of the radiative losses 
due to finite optical thickness and the increase due to the in- 
corporation of additional line emissions, would determine the 
actual shape of the loss function at low temperatures and, there- 
fore, the actual value of the growth rate of the unstable thermal 
modes. Also, the growth rates obtained here may be compared 
to the characteristic time scale for atomic processes. This means 
that the assumptions of ionization equilibrium of the plasma and 



electron thermal distribution might not be strictly valid at these 
short time scales. We plan to further investigate this aspect in the 
near future. 

We conclude that thermal instabilities may take place in 
prominences at lower temperatures than those predicted with 
previously existing loss functions. This may be important for the 
dynamics and energy balance of the prominence plasma. The ob- 
tained growth rates suggest that this low-temperature instability 
may have an observable effect in prominences. For example, this 
low-temperature instability may help to form density enhance- 
ments in regions where the plasma is already cool as in, e.g., the 
prominence threads observed in Ha images. However, nonlinear 
studies beyond the present normal mode approach are needed to 
assess the actual impact of the instability on the evolution of the 
prominence medium. This is relegated to future works. 

Acknowledgements. The authors acknowledge discussion within ISSI Team on 
Solar Prominence Formation and Equilibrium: New data, new models. RS ac- 
knowledges support from a Marie Curie Intra-European Fellowship within the 
European Commission 7th Framework Program (PIEF-GA-2010-274716). RS 
and JLB acknowledge the financial support received from the Spanish MICINN 
under project AYA201 1-22846. RS and JLB also thank the financial support 
from CAIB through the "Grups Competitius" scheme. SP acknowledges the 
support from the Belgian Federal Science Policy Office through the ESA- 
PRODEX programme. CHIANTI is a collaborative project involving George 
Mason University, the University of Michigan (USA), and the University of 
Cambridge (UK). 

References 

Barcelo, S., Carbonell, M., & Ballester, J. L. 2011, A&A, 525, A60 

Braginskii, S. I. 1965, Rev. Plasma Phys., 1, 205 

Carbonell, M., Oliver, R., & Ballester, J. L. 2004, A&A, 415, 739 

Carbonell, M., Terradas, J., Oliver, R., & Ballester, J. L. 2006, A&A, 460, 573 

Carbonell, M., Oliver, R., & Ballester, J. L. 2009, New A, 14, 277 

Cargill, P. J., & Hood, A. W. 1989, Sol. Phys., 124, 101 

Dahlburg, R. B. & Mariska, J. T. 1988, Sol. Phys., 117, 51 

Dere, K. P., Landi, E., Mason, H. E., Monsignori Fossi, B. C, Young, P. R. 1997, 

A&AS, 125, 149 
Field, G. B. 1965, ApJ, 142, 531 

Forteza, P., Oliver, R., Ballester, J. L. & Khodachenko, M. L. 2007, A&A, 461, 
731 

Gouttebroze, P., & Labrosse, N. 2009, A&A, 503, 663 
Heyvaerts, J. 1974, A&A, 34, 65 
Hildner, E. 1974, Sol. Phys., 35, 123 

Kaastra, J. S., Mewe, R., & Nieuwenhuijzen, H. 1996, in UV and X-ray 
Spectroscopy of Astrophysical and Laboratory Plasmas: Proceedings of the 
Eleventh Colloquium on UV and X-ray (Tokyo: Universal Academy Press), 
ed. K. Yamashita, & T. Watanabe, 41 1 



6 



R. Soler et al.: Thermal modes in prominences 



Karpen, J. T., Antiochos, S. K., Picone, J. M., Dahlburg, R. B. 1989, ApJ, 338, 
493 

Klimchuk, J. A., & Cargill, P. J. 2001, ApJ, 553, 440 

Labrosse, N., Heinzel, P., Vial, J.-C, Kucera, T., Parenti, S., Gunar, S., 
Schmieder, B., & Kilper, G. 2010, Space Sci. Rev., 151, 243 

Landi, E., Del Zanna, G, Young, P. R., Dere, K. P., & Mason, H E. 2012, ApJ, 
744, 99 

Lighthill, M. J. 1960, Phil. Trans. R. Soc. A, 252, 397 

Mazzotta, P., Mazzitelli, G, Colafrancesco, S., & Vittorio, N. 1998, A&AS, 133, 
403 

Milne, A. M., Priest, E. R., & Roberts, B. 1979, ApJ, 232, 304 
Oliver, R., & Ballester, J. L. 2002, Sol. Phys., 206, 45 
Oran, E. S., Mariska, J. T., & Boris, J. P. 1982, ApJ, 254, 349 
Parenti, S., Buchlin, E., Cargill, P. J., Galtier, S., & Vial, J.-C. 2006, ApJ, 651, 
1219 

Parenti, S., & Vial, J.-C. 2007, A&A, 469, 1109 
Parker, E. N. 1953, ApJ, 117, 431 

Priest, E. R. 1982, Solar Magnetohydrodynamics (D. Reidel Publishing 
Company) 

Rosner, R., Tucker, W. H. & Vaiana, G. S. 1978, ApJ, 220, 643 
Schure, K. M., Kosenko, D., Kaastra, J. S., Keppens, R., & Vink, J. 2010, A&A, 
508,751 

Soler, R. 2010, PhD Thesis, Universitat de les Hies Balears, Spain 
Soler, R., Oliver, R., & Ballester, J. L. 2008, ApJ, 684, 725 
Soler, R., Oliver, R., & Ballester, J. L. 2009a, ApJ, 699, 1553 
Soler, R., Oliver, R., & Ballester, J. L. 2009b, ApJ, 707, 662 
Soler, R., Oliver, R., & Ballester, J. L. 2010, A&A, 512, A28 
Soler, R., Ballester, J. L., & Goossens, M. 2011, ApJ, 731, 39 
Van der Linden, R. A. M., Goossens, M., & Goedbloed, J. P. 1991, Phys. Fluids 
B, 3, 866 

Van der Linden, R. A. M., & Goossens, M. 1991, Sol. Phys., 134, 247 
Zaqarashvili, T. V, Khodachenko, M. L., & Rucker, H. O. 201 1, A&A, 529, A82 



7 



